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Abstract 

A random phase property establishing a link between quasi-one-dimensional random 
Schrodinger operators and full random matrix theory is advocated. Briefly summarized 
it states that the random transfer matrices placed into a normal system of coordinates 
act on the isotropic frames and lead to a Markov process with a unique invariant measure 
which is of geometric nature. On the elliptic part of the transfer matrices, this measure is 
invariant under the full hermitian symplectic group of the universality class under study. 
While the random phase property can up to now only be proved in special models or in a 
restricted sense, we provide strong numerical evidence that it holds in the Anderson model 
of localization. A main outcome of the random phase property is a perturbative calculation 
of the Lyapunov exponents which shows that the Lyapunov spectrum is equidistant and 
that the localization lengths for large systems in the unitary, orthogonal and symplectic 
ensemble differ by a factor 2 each. In an Anderson-Ando model on a tubular geometry 
with magnetic field and spin-orbit coupling, the normal system of coordinates is calculated 
and this is used to derive explicit energy dependent formulas for the Lyapunov spectrum. 

PACS: 72.15.Rn, 73.23.-b, 73.20.Fz 
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1 Introduction 

The quantum mechanics of an electron in a disordered wire or a mesoscopic system is described 
by a random Hamiltonian which, in the tight-binding approximation, is a matrix typically given 
as the sum of a fixed kinetic part H and a random part \H\ containing a coupling constant 
A > 0. The location of the random and deterministic matrix entries then reflects the spacial 
structure of the sample. A widely used paradigmatic model is the Anderson Hamiltonian [And]. 
Calculating analytically physical quantities such as the density of states (DOS), the localization 
length, or the transmission amplitudes from a given random Hamiltonian is a difficult endeavor. 
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It usually becomes feasible only under the assumption that the randomness of the Hamiltonian 
somehow leads to one of the classical random matrix ensembles. The latter have invariance 
properties which allow to calculate the relevant averages. The most bold such assumption 
states that the Hamiltonian of the sample itself is drawn from a random matrix ensemble, i.e. 
the scattering matrix is given by one of the circular ensembles. This may be justified when 
studying universal properties both in the metallic phase and the localization phase but it is 
not valid in situations where spatial structure is important or where one wants to analyze a 
dependence on parameters of the model. 

In this work, we consider quasi-one- dimensional random systems described by random trans- 
fer matrices S = S(X) which may or may not be obtained from a random Hamiltonian Hq+XHi. 
The main aim is then to establish a natural link to invariant ensembles of random matrix theory. 
In a strictly one-dimensional case, this link degenerates to the random phase approximation 
[ATAF] which amounts to supposing that the phases of the Dyson-Schmidt variables |Dys[ [SchJ 
(called Priifer phases in the mathematical literature) are distributed according to the Lebesgue 
measure on the unit circle, which is also the unitary group U(l) of dimension 1. As is well-known 
[PFJ, this only holds in the regime of weak disorder and if one uses the correct coordinates, 
namely, the modified Priifer variables. Our non-commutative generalization to a situation with 
L channels uses matrix-valued Priifer phases taking values in a unitary group. These unitaries 
are in one-to-one correspondence with the isotropic frames (orthonormal coordinate systems 
within a Lagrangian plane w.r.t. the hermitian symplectic form) which hence carry a natural 
measure. This measure is, moreover, invariant under the natural action of the hermitian sym- 
plectic group on the isotropic frames (see Section [2] for details) . Now, the transfer matrices 
5(A) are in this hermitian symplectic group and certain cocycles associated to the action give 
the Lyapunov spectrum, a fact at the base of the transfer matrix method \MK\ IPS] which is the 
prime numerical procedure used for the calculation of the Lyapunov spectrum. Our random 
phase property (RPP) describes the distribution of the isotropic frames in the weak-coupling 
limit A — > when the transfer matrix is brought into a hermitian symplectic normal form: 

M^SM = TZe xv . (1) 

Here Ai is a symplectic basis change to the normal form TZ which is a rotation matrix on the 
elliptic channels and a simple linear expansion on the hyperbolic ones, and V is the random 
perturbation. Then the RPP states that the invariant measure on the isotropic frames in this 
coordinate system is given by the above invariant measure on the elliptic channels, combined 
with a deterministic distribution on the hyperbolic channels. The precise formulation of the 
RPP is given in Section [31 Let us note that one also finds the terms open and closed channels 
instead of elliptic and hyperbolic channels in the literature [Ben]. 

There is a different widely used non-commutative generalization of the one- dimensional ran- 
dom phase approximation going back to the works [Pol] IMPK] . It supposes that the unitary 
matrices in the polar decomposition of mesoscopic blocks of transfer matrices are distributed 
according to the Haar measure. This is referred to as local maximum entropy Ansatz (MEA) 
(the global one refers to the transfer matrices of the whole sample) or also isotropy assumption 
[Ben] . In the weak disorder regime it allows to deduce the DMPK equation for the flow of 
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transmission amplitudes [Dolt IMPK] . The latter lead to reasonable physical predictions of uni- 
versal nature, but without any parameter dependence (such as energy and coupling constants). 
There are also numerous physics papers discussing the validity of (variants of) the MEA itself 
(see |MS[ IMSt[ \MT\ \MC\ \CB\ IFYMS] and many others) , but we are unaware of numerical tests 
of its validity for concrete models. 

The basic but crucial difference between our RPP and the local MEA is the following 
(see Section H] for details): in the MEA the distribution of the transfer matrices themselves 
is supposed to be of maximal entropy, whereas in this paper the distribution of the transfer 
matrices is given and fixed by the model under study and then only the associated random 
dynamical system on the space of isotropic frames is supposed to have an invariant measure of 
maximal entropy (that is, a Haar measure). The RPP is considerably weaker than the MEA 
(see Section HJ), but it nevertheless allows to deduce the Lyapunov spectrum and its energy 
and model dependence. Of course, we recover the universal features (here the equidistance 
of the Lyapunov spectrum and the dependence of the localization length on the universality 
class). Moreover, the RPP provides a concrete statement that is easily falsifiable by numerical 
analysis. Numerics also allow to estimate the (mesoscopic) length scale beyond which the RPP 
is valid. We dare to call the RPP a property for the following reasons: The RPP was rigorously 
proved to hold in the Wegner L-orbital model |SS3j . Furthermore, for an Anderson model on 
a strip, which is precisely the model studied later in the introduction and in Sections [7] to 
[10] below, rigorous proof could be provided [SBlj that at least the low moments coincide with 
those calculated from the RPP (this is sufficient for the calculation of the Lyapunov exponents). 
Moreover, we will provide numerical evidence that the full RPP does hold for this model. 

The main application of the RPP presented here concerns the perturbative calculation of the 
L positive Lyapunov exponents 71 > . . . > 7^ > 0. Their definition is recalled in Section [2j The 
localization length is the inverse of the smallest positive Lyapunov exponent 7^. In Section [6] 
we prove a general perturbative formula for the Lyapunov exponents (Theorem [1]) whenever the 
RPP holds. It shows that the Lyapunov spectrum is equidistant and that the signature of the 
universality class (orthogonal, unitary or symplectic) are characteristic factors 2 between the 
localization lengths, in agreement with |Dol[ IEF[ IMCj . As these universality classes are closely 
linked to the complex, real and quaternion number field respectively [Mehj . we will consistently 
index various objects by a number field K which is either C, R or H. 

Rather than stating Theorem [TJ in this introduction, let us highlight the corollaries for the 
Anderson- Ando Hamiltonian [Andl lAndo] on a discretized tube. The Hilbert space describing 
a particle with spin on the tube is £ 2 (Z, C L ) Cg> C 2 where we think of C L as discretized annulus 
and C 2 is the spin degree of freedom. On £ 2 (Z, <C L ) <8> C 2 acts the right shift Si on Z and the 
cyclic shift S 2 on the fiber C L , namely {S2) L = 1. Neither of the shift operators S± and £2 
effects the spin degree of freedom. The magnetic spin-orbit Laplacian (in Landau gauge) is a 
bounded operator on £ 2 (Z, C L ) (g> C 2 given by 

H = S 1 + S* 1 +e^S 2 + e- t(p S; + 2tt{S 1 -S* l )®s y + 2tt{S 2 -S;)®s x (2) 

Here ip G [0, 2n) is half the magnetic flux through a cell of the lattice, t is a coupling constant of 
the spin-orbit coupling and the spin matrices s = (s x , s y , s z ) are as usual given by s x = ~ ( 5 )> 
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Furthermore, the random potential of Anderson type is given 
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(3) 



w n j \n,l)(n,l\ , 

neZ 1=1 

where the w n> i are independent and identically distributed centered real variables with unit 
variance. For sake of simplicity, we suppose the support of the distribution of the w n j to 
be compact. The Anderson- Ando Hamiltonian is then H = H + XHx- It depends on three 
parameters (p, t and A. For p ^ 0, it is in the unitary universality class with no time-reversal 
symmetry. If (p — and t — 0, the spin degree of freedom is effectively suppressed and H 
is in the orthogonal universality class of time-reversal invariant systems with no or even spin. 
Finally, if p = and t ^ 0, the model is in the symplectic universality class of time-reversal 
systems with odd spin. 

As usual \PF\ [SB2J the (not necessarily square integrable) solutions ip of the Schrodinger 
equation Hip = Eip at energy E can be calculated using transfer matrices. For the Anderson- 
Ando Hamiltonian and at height n e Z, these transfer matrices are 

[E - (S 2 + S%) - t(S 2 - S^)2is x - Xw n ](l + 2Us y )- 1 -(1 + 2tis v )* ' 

(i + 2ti S y)- 1 o 



s 
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where w n = diag(w nj i, . . . , w ni i)- This is a 2L x 2L matrix with 2x2 matrix entries describing 
the spin degree of freedom. In Sections [7] to [91 we will construct the hermitian symplectic 
basis change Ai needed to bring these transfer matrices into a normal form ([1]). Both the 
basis change Ai and the number L e of elliptic channels in the normal form TZ depend on the 
energy E, as do the Lyapunov exponents. Then, supposing that the RPP holds, one can use 
the general formula given in Theorem [1] in order to calculate the Lyapunov spectrum. The 
result is that 
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where ki = arccos(-| — cos(^)) for I = 1, . . . , L whenever the argument of the arcus sine is of 
absolute value less than 1. The number of Vs for which this is the case is precisely the number 
L e of elliptic channels and the sum in ([5]) carries precisely over these indices. A similar formula 
has already been derived by different means in [Do2j . For the inverse localization length 7*, 
this formula and its rigorous proof can already be found in our prior works fSBl} IRSj . and 
the case L = 1 has been known for a long time |PF| . Now for the case without time-reversal 
invariance, let us take t = and take the magnetic flux p to be small, but non- vanishing. Then 
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Finally, for the result in the symplectic universality class one takes tp = and t small: 
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Figure 1: The energy dependence of the smallest three Lyapunov exponents 7^, 7^_! and 7^_ 2 
for K = C (grey lines) and K = R (black lines) as given by ([HD and (jSj) respectively for a 
tube width L = 20 and disorder strength A = l.ll/\/l2. For comparison, numerical values 
calculated by the transfer matrix method (at 1% error corresponding to one standard deviation) 
are plotted for 7^ (o), 7 f 9 (□) and 7f 8 (o). 



These formulas can be compared to ab initio results obtained with the transfer matrix methods. 
Figure Q] shows that the agreement is remarkably good. Note that each of these formulas have 
singularities when ki is close to 0. This happens at so-called internal band-edges |RSj . The 
numerics of |RSj and Figure [1] show that these singularities are washed out. For an analysis of 
this phenomenon one needs to go beyond the RPP. 

The first important implication of formulas (JSJ), ([6]) and ([7]) is that the elliptic (bottom) 
part of the Lyapunov spectrum is equidistant with the same spacing. This equidistance also 
follows from the DMPK equations |MCj and is an important ingredient, for example, the 
analysis of shot noise [BBJ. On a more analytical level, Dorokhov [Do3j has already shown such 
equidistance for the Wegner L-orbital model (which is closer to random matrix theory because 
Hi is a full hermitian random matrix), a fact that was rigorously proved in |SS2j . The second 
important conclusion concerns the localization lengths, namely in the limit of large L e , we find 

\ll = 1 C L = 2jf. 

The first equality reflects the suppression of backscattering in presence of magnetic fields. These 
formulas agree with prior findings by different methods |EFl IDoll IMC] . 

As was already stressed in |RSj . these formulas hold only in a quasi-one-dimensional situa- 
tion where the sample is much longer than wide. All formulas have the same scaling predicted 
by Thouless |Thoj , namely as the sum over I is of order L e , one has roughly 7^ ~ UL-p+1). 
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There is basically only a difference of a factor of 2 between the localization lengths in the three 
universality classes, whereas it is well-known that the two-dimensional behavior is dramatically 
different (localization for R, derealization for HI and Landau levels for C) . This indicates that 
these two-dimensional behaviors cannot be analyzed perturbatively. Finally, let us direct the 
interested reader to Section [TU] where formulas for a quasi-one-dimensional Anderson model 
with a higher-dimensional fibers are derived. 

Before beginning the technical part of the paper, a short statement about the mathematical 
rigor may be at place. We definitely hope that Sections [2] and [5] as well as the appendices 
satisfy high standards in this respect, and also the diagonalizations in Section [7] to [9] are 
rigorous. Clearly the RPP itself is not proved here, and therefore the proofs of its implications 
such as Theorem [I] have only an algebraic value. However, using these algebraic identities the 
methods of |SBlj allow to prove ([5]) to (jTJ) rigorously up to an unknown density matrix (see 
[RSJ, with the RPP this matrix is the identity), albeit away from the internal band edges and 
with unsatisfactory error estimates. The main short-coming of the RPP (also on the level of 
theoretical physics) is the lack of a detailed analysis of the internal band edges |RSj . see the 
discussion in Section [3l Therefore the error estimates in (jSJ) to (0) and Theorem [T] break down 
at these internal band edges. 

Acknowledgements: This work was supported by the DFG. H.S.-B. thanks the Instituto 
de Matematicas in Cuernavaca for an extremely pleasant stay during his sabbatical. R.A.R. 
thanks A. Stuart and G. Golub for discussions. 

2 Action on isotropic frames and flags 

Let us first introduce the symplectic form J , the Lorentz form Q and the Cayley transform J 
as the following 2L x 2L matrices (matrices of this size are denoted by mathcal symbols in this 
work), each composed by 4 blocks of size L x L: 

The following identities will be useful: 

CJC* = ~Q, CJC* = - J. (8) 

i % 

Next let us recall that the quaternions H are the real span of four units q = 1, qi, q 2 and q 3 
which satisfy the algebraic relations q\ — q\ — q\ — q\qiq% = —qo- We will think of these units 
to be the following complex 2x2 matrices: 

90 = ( I 1 ) ' qi = ( -J ' 92 = ( -1 I ) ' q3 = ( * ) ' 

Hence a quaternion aGHisa2x2 matrix of the form a = Y^j=o a j c lj whh o,j £ K. Often 
q will be omitted. The quaternion conjugation is defined by a* = a — Yl^=i a j ( lj an d this 



6 



coincides with the adjoint of a as 2 x 2 matrix. A quaternion matrix A G Mat(L x L, H) is 
simply a matrix with quaternion entries, hence also a 2L x 2L matrix with complex entries. 
The adjoint A* of this complex 2L x 2L matrix is defined as usual. Note that A* is also given 
by the transpose of the quaternion conjugate of A. 

Even though simply expressed in terms of q 2 , it will be convenient to introduce I = ( ° ~q ). 
The reason is that we will allow / to act on arbitrary complex matrices and that it will be used 
to implement the time reversal symmetry below. Now a given a G Mat (2 x 2, C) is a quaternion 
if and only if a = Pal (here the overline denotes the complex conjugation). Further, somehow 
abusing notation we also write I = diag(J, ...,/) with as many entries as needed in a given 
situation. Hence I always has an even number of rows and columns and satisfies I 2 = — 1. Now 
given A G Mat(2L x 2L, C), one has the equivalence of A G Mat(L x L, H) and I*AI = A. 

Next let us define the hermitian symplectic groups HS(2L, K) for K = R, C, H: 

HS(2L,K) = {T G Mat(2L x 2L,K) \ T*JT = J } . 

In the case K = R the hermitian symplectic group HS(2L,R) coincides with the symplectic 
group SP(2L, R), but in the other two cases the use of the adjoint instead of transpose makes 
a difference. In our prior work |SB2j we nevertheless denoted these groups by SP(2L, K) (with, 
moreover, L replaced by ^ in the case K = H). As this leads to notational conflicts with the 
literature, we now plan to use the notation HS(2L, K). Let us note that 

HS(2L, R) = [T g HS(2L, C) \ T = T } , 

and similarly HS(2L, H) is a subgroup of HS(4L, C): 

HS(2L,H) = {Tg HS(4L,C) | I*TI = T } . 

The equations T = T and I*TI = T can respectively also be written as T l JT = J and 
T*/J7"T = I J, which were used as definition in fSB2j. Let us also note that HS(2L, H) is 
also isomorphic to the classical group SO*(4L). For our purposes below, it will at times be 
convenient to use the Cayley transform of these groups 

U(L,L,K) = CHS(2L,K)C* . 

These groups are called generalized Lorentz groups of signature (L,L), also called pseudo- 
unitary groups [MPK] . In fact, due to (jSD, U(L, L, C) is the group of matrices T G Mat(2L, C) 
satisfying T*QT = Q and hence conserving the Lorentz form Q. In U(L, L, R) they further 
satisfy either QJTQJ = T (or alternatively T*J7"T = J), while U(L, L, H) is the subgroup 
of U(2L,2L,C) characterized by PQJTQJI = T (or alternatively T l IJT = IJ). Unfortu- 
nately, U(L, L, R) and U(L, L, M) are not real and quaternion matrices other than the notations 
suggest. This can be considered a disadvantage of the groups U(L,L,K), on the other hand 
the associated frames are simpler than those in the hermitian symplectic representation. 

One can view J as a sesquilinear (hermitian) form on C 2L . A subspace of C 2L is called 
isotropic (w.r.t. J") if this sesquilinear form vanishes on it, that is v*Jw = for all vectors 
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v,w in the subspace. The maximal dimension of an isotropic subspace is L and such maxi- 
mal isotropic subspaces are also called hermitian Lagrangian. By definition, each element of 
HS(2L, C) maps an isotropic subspace to an isotropic subspace. A flag of isotropic subspaces 
is an increasing and maximal sequence of isotropic subspaces. The flag manifold F(L, C) is by 
definition the set of flags of isotropic subspaces w.r.t. J . Now F(L, M) is the subset of ¥(L, C) 
composed of flags satisfying v l Jw = for all vectors v, w in the flag. Similarly F(L, H) is the 
subset of F(2L, C) with flags satisfying v l IJw = for all vectors v,w. The group HS(2L,JK) 
naturally acts on F(L, IK). 

It will be very convenient to cover F(L, IK) by a simpler object. Each flag in F(L, K) can be 
described by a 2L x L matrix $ = (0i, . . . , (f>i) by letting the pth isotropic subspace be spanned 
by the vectors 0i, . . . , <p p . Moreover, one can choose these vectors to be orthonormalized. Hence 
let us introduce the set of (maximal) isotropic frames w.r.t. J by 

1{L, K) = {$ G Mat(2L x L, K) | $*$ = 1 , <F J<f> = } . (9) 

For sake of concreteness, let us note 

I(L,R) = {$ G I(L,C) |$ = $} , I(L,U) = {$Gl(2L,C)|r$/ = $} . 

Now, we will also work with the flag manifolds C F(L, K) of isotropic flags w.r.t. the sesquilin- 
ear form Q satisfying possibly the real or quaternion symmetry. In the same way, CI(L,IK) 
contains the isotropic frames w.r.t. Q. These sets are particularly simple. In fact, because 
CI(L,C) are those L- dimensional frames in C 2L satisfying = 0, one readily verifies that 



CI(L,C) 



J_(U 



U,Ve\J(L)\. (10) 



This also allows to write out a formula for elements of I(L, C) in terms of two unitaries. 
Furthermore, the relations C*$ = C*$ and J*C*$/ = C*$ respectively lead to 



CI{L, 



-(- 



U G U(L) } , CI(L, 



U G U(2L) 



(11) 

Because of these representations there are natural measures on C 1(L, K) (and hence also 
I(L, K)) induced by the Haar measures on the unitary groups. Let us note that under the 
stereographic projection $ = (£) G I(L, K) i— > (a — «6)(a + z6) _1 [SB2J these measures give those 
of Dyson's circular ensembles CUE, COE and CSE respectively on the maximally isotropic sub- 
spaces. We will next analyze the projection from I(L, IK) to F(L, K) and then these measures 
also lead to natural measures on F(L, IK). 

Two frames $, $' G I(L, C) describe the same flag if and only if $ = <&'S for un upper 
triangular matrix S. But $*$ = 1 and ($')*$' = 1 imply S*S = 1. The only upper triangular 
unitaries are the diagonal unitaries, which in turn can be identified with the torus T L . Hence 
we conclude that the isotropic frames form a T L -cover of the flag manifold: 

F(L,C) = I(L,C) /T L . (12) 
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For IK = M, the diagonal unitary S has to be real so that S G (Zy L , while for IK = H it must 
satisfy I*SI = S and the set of such diagonal unitaries is in bijection with the torus T L . Thus 

F(L,R) = I(L,R)/(Z 2 ) L , F(L,H) = I(L,H)/T L . 

Clearly the group HS(2L, K) naturally acts on F(L, K). Besides the simple expressions in 
( TTUT) and ( TTT1) . the main reason for introducing I(L, IK) is that the action of HS(2L,IK) lifts 
in a natural way which is particularly useful for the calculation of Lyapunov exponents. Let 
us first study the action of the unitaries in HS(2L, IK), or equivalently of the unitaries in 
U(L,L,K) on CI(L,K). For K = C, one has U(L, L, C) R U(2L) = diag(U(L), U(L)) and the 
action on CI(L,C) is simply given by left multiplication. By (TlO|) the measure on CI(L, C) 
induced by the Haar measure on U(L) x U(L) is invariant under this action. Furthermore, each 
T G U(L, L, R) n U(2L) = U(L) is of the form T = diag(W, F)_with W G U(L) and, similarly, 
T G U(L, L, H) nU(4L) = U(2L) is of the form T = diag(PF, JWJ) with W G U(2L) (this can 
be checked using Lemma 2 of |SB2J). Again, it follows from fTlT|) that the invariant measures on 
C I(L, IR) and C I(L, M) are invariant under these actions. Clearly this transposes to the actions 
of the unitaries in HS(2L,K) acting on I(L, IK). 

Now we want to extend the action of the unitaries in HS(2L,K) to a group action of all 
of HS(2L,K) on I(L,K). Given T G HS(2L,K) and $ G I(L,K), the 2L x L matrix T$ is 
isotropic w.r.t. J but not orthonormalized. However, applying a Gram-Schmidt procedure one 
then obtains a new isotropic frame T • $ G I(L, IK). In order to write this out more explicitly, 
let $ = (0i, . . . , 4>l) with (pi G IK 2i . Then T$ is isotropic w.r.t. J and the action is defined by 
$' = T ■ $ where = ((/>[, . . . , 0^) is calculated iteratively by 

^ = ' = T ^ p ~ K)* T 0p • ( 13 ) 

ll^pll ?=1 

Note that (0')*T0 P G IK and that the order of the factors in the last term is important in the 
case K = EL Now for T G U(L, L, K) and $ G CI(L, K), we define T ■ $ = C(C*TC • C*$). 

Proposition 1 Let T G HS(2L,K) and $ G I(L,K). 

(i) There exists a unique upper triangular matrix S(T, $) G Mat(L x L, IK) with positive entries 
on the diagonal such that 

r ■ $ = t^s^t,®)- 1 . 

(ii) The map (T, $) i— > S^T, $) a multiplicative cocycle for the action o/HS(2L, K) on 
I(L,K), 

namely if also T G HS(2L, IK) 

S(TT', $) = ^(T, T' • $) ^T', $) . 

(iii) T/izs cocycle does not project to a cocycle on F(L, IK), but one has for T G T L 

S(T,$T) = T- l S(T,$)T. (14) 
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Proof. Parts (i) and (ii) are immediate from the considerations above and for (1131) . one has to 
analyze (1T3|) a bit more closely. □ 

Let e p , p = 1, . . . , L, be the standard (real) basis of K L . Then e*S(T, $)e p > is the pth 
entry of on the diagonal of S(T, $) and one can define the additive real-valued cocycles 

g p (T,$) = r log( e ;5(T,$)e p ) , 

where r = 1 in the cases K = C, R, and r = | Tr2 in the case K = H if the quaternions are 
identified with 2x2 matrices (hence r just extracts the real coefficient of go)- It follows from 
(TL4|) that this is actually a cocycle on the flag manifold F(L, K). Note that in the case K = H, 
e*S(T, $)e p is a quaternion 2x2 matrix, but by Proposition CD the diagonal is real and hence 
each diagonal entry appears at least twice. It can be convenient to extract the diagonal entry 
by taking half of the trace of the quaternion. 

Now given a random process {T n ) n >i in HS(2L, K), one can define the growth exponents of 
these cocycles as usual by 

1 1 - 

7 ? = i im K9p(TN'--T lt $ ) = Mm - J2 9p{T n ,®n-i), (15) 

N^oo iV iV— >oo iV ' 

n=l 

where 

$ n = (T n ---T 1 )-$ = T n $ eI(L,K), (16) 

and $o is an initial condition. These growth exponents are then actually the usual Lyapunov 
exponents, namely one can check [SBlJ that for almost every initial condition $ = (</>!, . . . , <Pl) 

V7J = lim — log det p ((4>i\T N ---T l \4> k ) 1 <i }k < p ) . 

9=1 

Let us point out again that in the case K = H, (|T5|) only defines the L possibly distinct 
non-negative Lyapunov exponents. The Lyapunov spectrum is indeed twice degenerate and 
obtained by doubling each of these exponents [SS3J. 

3 The random phase property 

In many models of solid state physics one is naturally led to study the random action of 
hermitian symplectic transfer matrices on flag manifolds and the isotropic frames. Examples 
are provided in Sections [7] to [10J This leads to a Markov process on the isotropic frames of 
the type given in (1T61) . In a quasi-one- dimensional situation there are associated invariant 
measures on the isotropic frames and it is typically very difficult to determine them explicitly 
in an analytic manner. However, in a perturbative situation of weak and isotropic coupling of 
the randomness we believe that there is a unique such measure which is, moreover, close to 
a geometric invariant measure, provided one chooses normal coordinates. The random phase 
property (RPP) makes this claim precise and verifiable. Below we state in which cases and to 
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what extend the RPP can actually be proved and below we provide numerical evidence that 
the RPP holds in a variety of particular models. In Section [6] we show what kind of conclusions 
can be drawn from the RPP. 

We consider the following set-up. Let be given a random process ( / 7^) n >i hi HS(2L, K) in 
the normal form 

T n = TZe xv " , (17) 

where A > is a coupling constant and V n G hs(2L, K) are independent and identically dis- 
tributed Lie algebra elements, and TZ G HS(2L, IK) is a hermitian symplectic matrix of the form 
TZ = lZ e lZh = TZ}JZ e where TZh is an expansion matrix and TZ e G 0(2L) is a rotation matrix. In 
order to describe these matrices more explicitly, let k = diag(/«i, . . . , k l ) with diagonal entries 
Hi G K ordered according to their modulus and such that there is an Lh with I^lJ > 1 and 
Ki = 1 for / = L h + 1, . . . , L. Further be given a real matrix 77 = diag(r/i, . . . , r] L ) such that 
rji = for I = 1, . . . , Lh. Then the commuting matrices TZh and TZ e are given by 

_ / k \ _ / cos(r?) -sin(7/) 

We call Lh the number of hyperbolic channels, and L e = L — Lh the number of elliptic channels. 
Then let n e be L x L projection matrix of rank L e such that K,ir e = n e and let iih = 1l — TTe- 
One has L h = dim(7r /i ) and L e = dim(7r e ). We also set IT^ = diag(7T/ l , ir h ) and Il e = diag(7r e , n e ). 
They are called the projections on the hyperbolic and elliptic channels respectively. One has 
TZ e Yih = lift and TZhH e = n e . Here the qth channel, q = 1, . . . , L, is the span of e q , ei+ q G K 2L 
(where e g , q = 1, . . . , L or q = 1, . . . , 2L, denotes the standard basis of K L or K 2i ). Thus the 
qth. channel is elliptic or hyperbolic pending on whether 7i e e g = e q or n h e q = e g . 

In our examples in Sections [7| to [TU] below, quite some algebraic work is needed in order to 
bring the transfer matrices into the normal form (1171) . Furthermore, this normal form cannot be 
obtained for all energies, but at so-called internal band edges it would have to contain Jordan 
blocks. The use of the normal form is explained shortly. 

Now we consider the random action ffl6|) of the process (T n ) n >i on isotropic frames. Then 
($n)n>i is a Markov process, strictly speaking a family of Markov processes indexed by the 
coupling constant A. According to ffTU]) . there are unitaries U n and V n such that 

C®n = A=(Z n ). (18) 



In the cases K = R and K = H, one has V n = U n and V n = I*U n I respectively. The RPP 
describes the distribution of the $ n (or equivalently U n and V n ) based on the following heuris- 
tics. The hyperbolic rotation TZh dominates the random dynamics and forces the first frame 
vectors deterministically into the hyperbolic channels. The elliptic rotation also dominates 
the randomness, but only generates tori in I(L, K). The random perturbations V n change the 
orientation of the axis of these tori in a diffusive manner. The RPP now states that in a per- 
turbative limit first of all the splitting of hyperbolic and elliptic channels holds, second of all 
the hyperbolic frame vectors are ordered deterministically according the size of the hyperbolic 
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expansion factors and third of all that the random perturbation is maximally efficient in the 
sense that the distribution on the elliptic channels has maximal entropy and is given by the 
Haar measure. Thus the RPP states that the invariant measure of the Markov process $ n is 
unique and of geometric nature. Averages w.r.t. this invariant measure will be denoted E. In 
particular, it allows to calculate Birkhoff averages: 

1 - 

E /(*) = lim - /( $ «) . / e C(I(L,K)) . 

n=l 

In order to describe this invariant measure of the RPP, we consider Tr e Uir e , TiJJiih, ^eV^e etc. 
as matrices of size L e x L e , L e x L^, L e x L e , etc. 

Random Phase Property. In the limit A —>■ of weak coupling, the invariant measure of 
the Markov process on I(L, IK) generated by the normal form matrices (fTTj) is unique and, with 
errors of order 0(X 2 ), given by the following properties: 

(Rl) U = 7r e U7r e + n h Un h 

(R2) HhUiTh is deterministic and given by ir^ unless the moduli of tzi, I = 1, . . . , Lh, are degen- 
erate. 

(R3) Tx e Uix e is distributed according to the Haar measure on U(L e ) in the cases K = R, C, and 
to U(2L e ) in the cases K = H. 

(R4) In the case K = C, U and V are independent and identically distributed. 

Let us stress that by (Rl) 7i e Unh and 7thUn e vanish with errors of order C(A 2 ). Both (Rl) 
and (R2) can be proved by the techniques of |SBll Proposition 3] and [SB21 Lemma 7], with 
error estimates which are not optimal though. Item (R3) is the central piece of the random 
phase property. It establishes a connection with random matrix theory (here the Dyson circular 
ensemble). In the perturbative calculation of the Lyapunov exponents in Section [6] below, one 
actually does not need the full strength of the (R3), but only uses it to evaluate second and 
fourth moments of the unitaries as given in Appendix A. These identities for the second and 
fourth moments hold for a much wider class of distributions than those stated in (R3). The 
point is though that we believe (R3) to hold and give numerical evidence in Section [71 

A serious short-coming of the RPP as stated is that it does not cover the (vicinity of) 
internal band edges |RSj . Such internal band edges correspond to normal forms having Jordan 
blocks. When considering families of matrices (typically indexed by the energy), these Jordan 
blocks appear at points where an elliptic channel passes through a parabolic one to become 
hyperbolic channels (or inversely). When the RPP is applied to the calculation of the Lyapunov 
spectrum, this leads to singularities (see formulas <^ to ^7J) which are effectively smoothed out 
as show the numerics in |RSj . An rigorous analysis of this smoothing in the one-dimensional 
case is carried out in |SS1] . 
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4 Comparison of the RPP with the MEA 



The maximal entropy Ansatz (MEA) [DoTj IMPK| IMS] is a claim about the matrix entries of 
the following polar decomposition of the transfer matrix of a sample of length N (which in the 
local MEA is a building block of a longer sample): 

rn - i-r- f UN \ / y/1 + A N y/A^ \ ( U' N \ , . 

CTn "-^ C ~ \ v N ){ VT+An- ){ v' N ) ■ (19) 

Here A^v > is a diagonal L x L matrix with non-decreasing diagonal entries and u^, v^, 
u' N and v' N are in U(L). If K = R one has vn = and v' N = u' N while for K = R one has 
vn = I*unI and v' N = I*u' N I. The matrix is unique, but the polar decomposition is not 
because diagonal unitary factors (in T L ) may be added (and further unitary factors if A N is 
degenerate). The nicest feature of the polar decomposition is that the scattering matrix of the 
sample can be directly read off. 

The MEA now consists in supposing that for N sufficiently large the unitaries un, i>jv, u' n 
and v' N are Haar distributed and independent, apart from the correlation in the cases K = R 
and K = H described above. This corresponds to a system with only elliptic channels (in fact, 
a detailed treatment of closed channels is not known to us, but could easily be given along the 
lines of the last section). Hence (Rl) and (R2) become irrelevant and we now argue that the 
MEA implies (R3) and (R4) of the RPP. Indeed, using (113)) . 

1 ( Un \ r ^ _(u N 0\J_fu N \ 1 ( u N u N \ 
V2 V V N J N ~ V v n J ' V2 V v N J " 72 V V NV N ) ' 

where 

J_ f u N \ / VI +Aiv v 7 ^ \( u' N \ 
V2 V vn J " V v 7 ^ v/TTA^ J \ < ; ' ° ■ 
Thus, if and vn are Haar distributed (as they are by the MEA), so are Un and Vjv irrespective 
of what the distribution of mat and vn is. 

In conclusion, the MEA implies the RPP. In both approaches one makes an assumption 
on the invariant measure of a Markov process associated to a random dynamics generated by 
the transfer matrices. In the RPP the state space of the Markov process is compact and the 
dynamics relatively simple as described in Section [21 on the other hand in the MEA the state 
space is non-compact (namely HS(2L, K)) and the dynamics in the polar representation is com- 
plicated, that is, it is relative intricate but possible to calculate («/v+i, ^jv+i, Aat+i, u' N+1 , v' N+1 ) 
from T N+1 and (u N ,v N ,A N ,u' N ,v' N ). 

5 Perturbative formula for the additive cocycles 

Section [2] shows that the ergodic limits of the cocycles g p are precisely the Lyapunov exponents. 
As our main aim in Section [6] will be to develop a controlled perturbation theory for the 
Lyapunov exponents, the following proposition will prove to be helpful. It is actually quite 
easy to derive the lowest order terms of the expansion (and this was done in [SBlJ). However, 
the error estimate below is much better than previous ones. 



13 



Proposition 2 For V G hs(2L, C) and $ G I(L,C), we introduce the self-adjoint L x L 
matrices 

P x ($) = I <f>*(V + V*)<$> , P 2 ($) = i$*(2P*P + V 2 + (P*) 2 )$ . 
T/ien /or A G K. and the notation p p = YT q =i e q e q> 

a p (e AP ,$) = Are^^ep + AV [e;P 2 ($) ep +( e ;P 1 ($)e p ) 2 -2 e ;P 1 ($)p p P 1 ($)e p ] + 0(A 3 ) , 

with an error bound depending only on the norm ofV (in particular, independent of L). 

Proof. Let us begin by calculating upper triangular matrices S\,S 2 with real entries on the 
diagonal such that S(e xv , $) = 1 + XS ± + \ 2 S 2 + 0(\ 3 ). The definition of Pi($) and P 2 ($) 
is such that |e AP $| 2 = 1 + 2\P 1 ($) + 2A 2 P 2 ($) + 0(\ 3 ). Hence the defining property of 
S(e xv , $) leads to the equation 

1 = (l + XS^ + ^S;)' 1 [l + 2\P 1 ($) + 2\ 2 P 2 ($)](l + \S 1 + \ 2 S 2 y 1 + 0(\ 3 ) 
= 1 + A [2P 1 ($)-S 1 -S*} 

+ A 2 {2 P 2 ($) + S* 1 S 1 + - 2 P 1 ($)]S' 1 + - 2 P x ($)] - S 2 - S* 2 } + 0(\ 3 ) . 

Now we use the following basic fact: the unique upper triangular matrix S with real diagonal 
satisfying S + S* = P for some self-adjoint P is given by the strict upper triangle of P plus half 
the diagonal of P. Let S denote the corresponding super-operator, namely S = S[P]. With 
this notation, 

^ = S[2P!($)] , S 2 = ~{2P 2 ($) + SIS, + [S, -2P 1 ($)]S 1 + S^S? - 2P 1 ($)}} . 

For any triangular S and any analytic function /, one has e*f(S)e p = f(e*Se p ). Applying this 
to the logarithm and expanding gives 

g p (e xv ,<S>) = re;^XS 1 + \ 2 S 2 -^X 2 (S 1 ) 2 ^e p + 0(X 3 ) . 
Using e* S[P]e p = \ e*Pe p , some algebra now leads to 

g p (e xv , $) = At e^*)^ + y r {^^6,,+ [ e ;P 1 ($)e p ] 2 -4 Ke^Px^^eJ } , 

again with an error of order C(A 3 ). Finally, replacing the identities 

e;SlS iep = 4e;P 1 ($)p p _ 1 P 1 ($) ep +[ e ;P 1 ($)e p ] 2 = 4e;P 1 ($)p p P 1 ($)e p -3[e;P 1 ($)e p ] 2 , 
and 

e;P 1 ($)S 1 e p = 2e;P 1 ($)7r p P 1 ($) ep -[e;P 1 ($)e p ] 2 = Me[e;P 1 ($)S 1 e p \ , 
completes the proof. □ 
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6 Perturbative formulas for the Lyapunov exponents 



As an implication of the RPP, we show in this section how it allows to derive perturbative 
formulas for the Lyapunov exponents associated to a random process of the form ( flTj) . For 
sake of simplicity, we also assume a number of further properties. They are not essential, but 
simplify the algebra. Hence let us suppose that the random perturbation V is centered and 
that |/ci| > |/C2| > . . . > |klJ > 1 = >^L h +i = • • • = kl- Such strict inequalities hold generically. 
Then part (R2) asserts that ■K^U'Rh = tt^- 

The starting point is the formula (fT5l) for the Lyapunov exponents, but we suppose that the 
Birkhoff sum on the r.h.s. of (115j) can be calculated using the RPP which provides a measure 
on I(L, K). Denote the corresponding average combined with that over T by E, we therefore 
have 7^ = E<7 P (T, $). Appealing to Proposition [Qn), we deduce 

7p K = E g p (K h , e xv ■ $) + E g p (e xv , $) , (20) 

where we used that for the unitary lZ e one has g p (JZ e , $') = 0. The first term can be of order 
(9(1), while the second one is always of order 0(X 2 ) due to Proposition [2] and because V is 
centered. 

Let us begin by considering the first contribution in (120]) . For this purpose we need to 
know the upper triangular matrix S(lZh,®) where $ = | (^-v)) e I(L,K) is given in the 
representation (|TU|) and characterized by 

Now by (Rl) one has n e U = 7t e Uir e G U(L e ), etc., and 7Tft,C/ = n h and 717^ = n h , the latter also 
in the case K = R, H because n h is real and commutes with /. Thus S(TZh, $) = ir e S(lZh, $)vr e + 
ir h S(TZ h , ®)n h = 7T e + 7T^|fi;| = |«| and 

7 J = In(|«p|) + 0(A 2 ) . (21) 

For a hyperbolic channel p, this gives the lowest order term, while for an elliptic channel this 
lowest order contribution vanishes. 

We now consider the more interesting case of an elliptic channel p. It can be argued [SBlJ 
that the first term in (|20|) does not contribute even to order C(A 2 ), and therefore we now focus 
on the second term in (1201) . For this let us use the expansion formula in Proposition [2j Because 
V is centered, we have with errors of order 0(X 3 ), 

\2 / P \ 

7 J = — Er e* p <l>*(2V*V+V 2 + (V*) 2 )§e p - ^(2 - 5 q>p ) e* p $*(V*+V)<5>e q e* q $*(V*+V)<5>e p . 

(22) 

Here the expectation E is over the random Lie algebra element V and the distribution of the 
isotropic frames $. The main result of this section is that the Lyapunov spectrum given by 
this formula is equidistant at small coupling whenever the distribution of the $ is given by the 
RPP. 
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Theorem 1 Let the Lyapunov exponents for the L e elliptic channels be given by (1221) . that is 
p > Lh- Suppose the RPP holds. Then, with errors of order 0(X 3 ), 



X 2 



4L e (L e + 4 



2 4 



L-p + ^5 K , c + 8 KM + ^5 KM j ETr[LI, TP ! + V)\IVU, 



where 5 a ^ is the Kronecker delta equal to 1 if a = b and equal to otherwise and the trace Tr 
contains a factor | in the case K = H (similar as the | contained in r). 

According to (|22l) one has to calculate averages of certain functions on the isotropic frames. 
Instead of the self-adjoint operators 2V*V + V 2 + (V*) 2 and V* + V, let us consider for any 
self-adjoint A,B& Mat(2L x 2L, K) and indices 1 < p, q < L the quantities 



-p ) 



where the index K indicates the symmetry class. Let us also set A e = n e AfI e and B e = U e BU e . 
As J commutes with Il e , it follows B e G hs(2L, K) whenever B G hs(2L, IK). 

Lemma 1 Let A,B& Mat(2L x 2L, IK) be self-adjoint and suppose that B G hs(2L, K) and B 
is also in the Lie algebra of the special linear group. Let the RPP hold and let L h < p,q < L 
be elliptic. Then, up to errors 0(X 2 ): 



(iii) lf„{B) 



2 L P 



q'=l 



+ 



1+6.. 



p.q 



\L\ ^ AL e (L e + 1) 
Proof, (i) First let us write out $ in terms of U and V 



+ 



6 



p.q 



ALJ2L, - 1 



Tr[B 2 e 



1 



ETr 



u + v 

i{U -V) 



A 



U + V 
i{U-V) ) 6pe p 



Note that for K = M the trace contains a factor | stemming from the \ in the definition of r. 
Now because p is elliptic, by (Rl) we can replace A by A e with errors of order 0(\ 2 ). Thus 



A If (A) = ETr 



U* 



UcpCp 



ETr 



V* 



Ve p e* 



+ ETr 



U* 



Ap 



Ve p e* 



+ ETr 



V* 



Ap 



UepCp 



The first summand can now directly be calculated by Lemma [2] stated in Appendix A, first for 
the cases K = C, R: 



ETr 



U* 



-Tr 



Ap 



Tr [ e P <] 



Tr[A e (l+iJ)) 
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The second summand gives the same result except for the inverse sign. The sum of the first and 
second term thus give (i). Indeed, the third and fourth summand vanish, for IK = C because of 
the independence of U and V, and for K = M. because either the factor U appears twice or the 
factor U does, so that again the average over the unitary group vanishes. 

In the case K = H, the unitary group has dimension 2L e , but Tr(e p e*) = 2, so the end result 
is the same because Tr (I* AI) = Tr(A) for a quaternion matrix A. 

(ii) First let us note that Y^Li e q/ e* ql = ir h . Thus by (Rl) and (R2), $7r fe $* = ll[ h (l + jg), 
it follows immediately from (i) that (ii) = ^-Tt {U e BU h {l + J g)BU e ). But £> = £>*€ hs(2L, K) 
implies JB = BJ*. Therefore JQJ* = —Q combined with the cyclicity of the trace implies 
Tt(U e BU h jgBU e ) = -Tr(U e BU h jgBU e ) = so that the result follows. 

(iii) We now need to calculate fourth moments. Let us first expand in terms of U and V 
such that Ip q (B) can be written as 



— ETr 

16 



U + V 
i{U - V) 



B P 



u- 

i{U 



V 

-v) 



e p e p 



U + V 
i(U-V) 



B P 



U- 



V 

■v) 



e q e g 



Factoring U and V in each of the four factors, one has 16 terms. However, all but 6 of them 
vanish because one needs the arguments to come in complex conjugate pairs (U, U) and (V, V) 
in order to have a no n- vanishing average over the unitary group (cf. Appendix A). Even though 
very lengthy, let us write out the result: 
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ETr 



U* 



+ ETr 



+ ETr 



+ ETr 



+ ETr 



+ ETr 



V* 



v* 



u* 



v* 



B„ 



1 

—i 

1 
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1 

—i 

1 

i 

1 

—i 



B P 



B P 



B P 



B, 



B P 



Ue p e* p U* 



B P 



1 

—i 



Ve p e* p V* 



1 

—i 

1 

i 



Ue p e* p V* 



Ve p e* p U* 



Ve p e* p V* 



Ue p e;U* 



i 

1 

i 

1 

—i 

1 

i 



B„ 



B P 



B P 



B P 



Ue q e* 



B P 



Ve q e q 



Ve q e* q 



Ue q e q 



Ue q e* q 



Ve q e* q 



(23) 
(24) 
(25) 
(26) 
(27) 
(28) 



Now each of these terms can be calculated using Lemma [2J The matrices B and D in Lemma [2] 
are always given by B = B l = e p e* and D = D l = e q e q , while A and C differ in each of 
the terms. If we assume that p ^ q, one has BD = BD l = so that many terms vanish. 
Furthermore Tr(B) = Tr(D) = 1 if K = C, R and Tr(B) = Tr(7J) = 2 if K = H. We will show 
that (T23]) to (J26]) vanish, and that fT27|) = (T281 give the contribution on the r.h.s. in all cases. 

We now first focus on the case K = C and p ^ q. We deal with fT23l) using Lemma EJ^iii): 



LI 



-Tr 



B P 



B„ 



L e (L 



-Tr 



B P 
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Hence 



-^-^r [B\ - B e JB e j] 



-Tr[B e (l+iJ)) 



L e {Ll - 1) 

However, the last term vanishes because Tr(£> e ) = Tr(J'Be) = 0, and JB e J = B e J*J = B e 
implies Tr[B e J'B e J'} = Tr[i3g] so that also the first summand vanishes. As V has the same 
distribution as U, (I24p gives the same result as (I23p except for the minus sign in front of the 
imaginary units and it thus follows (124"]) = 0. Furthermore, (125]) is calculated by iterating 
Lemma [2]^i) twice. Thus 



d25D 



Tr 



ETr 



V* 



1 

—i 



B e 



1 

—i 



T T * * 

v e q e q e p e p 



Thus for p ^ q, one immediately sees (125]) = 0, but actually this always holds because the 
first factor vanishes for the same reason as above. Similarly, (126]) = 0. Finally, again using 
Lemma [2]^i) twice, 



(ETD 



LI 



Tr 



1 

—% 



B e 



B P 



Tr[£>g + B e J B e J] 2Tr[£ e 2 



LI 



As ( 127]) = (128]) . adding up all terms leads to the prefactor for K = C and p ^ q. 

For p = q each of the six terms above has two supplementary contributions by Lemma [2j 
Let us start with: 



L e (L e + 1) 



Tr 



B e 



n 2 



Tr 



Bp. 



B F . 



As each of the two traces vanishes by the calculation above, (T53D = and similarly fliM]) = 0. 
We already showed (l2"5l) = (]2"6"]) = 0. Finally, (TSTj) and f]2"g]) are the same as in case p ^ q and 
thus (iii) for IK = C is shown also for p = q. 

Now K = R, which allows us to use for the real and self-adjoint B the identity 



B P 



B P 



(29) 



The evaluation of (123]) and (1!M]) is the same as in the case K = C and thus (123]) = (124]) = 0. 
Furthermore Lemma [2]^iv) implies that (1251) = (I26p = because BD = BD l = for p ^ q, 
while for p = q 

<m = m = T (T l , 7\ (Tr [B 2 e - BeJBeJ] + Tr \B e f + Tr \JB e f) , 
which vanishes again. For ( 1271) and ( 128]) we will use Lemma E](v): 



d2TD = d2HD 



L e {L e + 1) 



Tr 



B P 



Be 



2Tr[i3g] 
^ e (^ e + 1) 
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which gives the case K = R for p ^ q. For p = q, there are two more contributions in 
Lemma [2^v) showing that both (1271) and (|27|l are exactly twice as large, showing the formula 
also for p = q. 

In the case K = H, the changes w.r.t. the case K = R are that now V = I*UI and 
Tr(e p e*) = 2, and that 2L e replaces L e as the dimensions on the unitary group. Moreover, we 
have at our disposal that I*BI = 1*3*1 = B for the self-adjoint quaternion matrix B so that 
( |29l) holds with factors of / and I* from left and right on either one of the identities. The 
verification of (1251) = ff2M = is the same as in the case R, that of = = also modulo 
the use of the modification of identity (|29|) . In the evaluation of (I2T1) and (I28I) . all four terms 
of Lemma E](v) contribute. With some care one gets the corresponding prefactor. □ 

Remark. As a check of the above formulas it is worth noting that the sum rule 



E = tt 1 ' (30) 

q'=L h +l e 

can be checked independently starting from the identity 

V $e,e*,<T = $tt = -n + ic* f u ^v*\ c 



The contribution corresponding to |ll e is evaluated using Lemma [T](i) and gives the r.h.s. of 
( 130]) . because the contribution of the second summand vanishes (as shows some further algebra). 
Then one has, in all three cases of Lemma [T](iii) indeed (as it should be by Lemma [1]) 

Tl '^e] _ (T _ -H rK , tK / 

4 ^ ~~ ^ e J p,<? + p,p ' vt q ■ 

Proof of Theorem [TJ, We need to evaluate each term in (1221) . The first term is calculated 
using Lemma [T](i) and the identity V* = JVJ: 

lf[2V*V + V 2 + (V*) 2 ] = ^Ti[U e (V + V*)U e VU e ] + ^ Ti[U e (V + V*)U h VU e ] . (31) 

Lie L>e 

The sum over q is split into a hyperbolic part 1 < q < and an elliptic one + 1 < q < L. 
The hyperbolic one is given by 



Y, 2I f, q ( V * + V ) = T Tr[U e (V + V*)U h VU e 



9=1 

where the equality follows from Lemma [H(ii) and again V* = JVJ. Due to the sign in (122|) . 
this hyperbolic part thus cancels with one of the terms on the r.h.s. of (l3~Tj) . The elliptic part 
is given by 

p 

E ( 2 - <W J p K ^* + P ) = 2{ P - ^ - !) J p K p'(^ + ?) ~ 2 J p K p(^* + V) » 
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where p' ^ p. Now only remains to evaluate the contributions using Lemma [TJiii) and then 
carefully resemble all terms. Let us exemplify with the case K = R. 



IfJV + V) = 2 



so that 



2 • 2 



[L e 4L e (L e + 1) 



(p-L h -l) 



1 2-2-2 

2 4L e (L e + l)J 



which leads to the formula in the case IK 



□ 



7 Magnetic Anderson model for tube geometry 

In order to calculate the Lyapunov spectrum more explicitly than in Theorem [TJ one needs 
to include further model specific information. From this section on, we will study the one- 
particle tight-binding Hamiltonians H = H + Ai^i described in equations (j2J) and ([3]) in the 
introduction. In this section, we take if ^ and t = 0. Hence the spin degree of freedom is 
suppressed and can thus be neglected. Therefore the transfer matrices (j3J) at energy E e R 
reduces to an 2L x 2L matrix 



S 



El- (e^S 2 + e'^SZ + X w n ) -1 
1 ' 



(32) 



where the dependence on n is suppressed. The transfer matrix S is in the hermitian symplectic 
group HS(2L, C). We first bring the transfer matrix into a normal form. Let us introduce, for 
Z = 1,...,L, 



Si 



( Kl \ 

\ fl,L ) 



e C 1 



fl,k 



1 



cxp 



Itxi Ik 
L 



m 



(./: 



i- 



/x,). Hence m is the discrete Fourier transform (in U(L)) and 



Then we set 

= diag(m,m) e HS(2L,C) n U(2L). One has m*^ 2 m = e ir > with = ^diag(l, 2, . . . ,L). 



Therefore m* (E - e ltp S 2 - e~ ltp S* 2 



E — 2 cos(ip + rj) is diagonal as well and is non-degenerate 



(for <p ^ 0). Let q be an adequate L x L permutation matrix such that the diagonal entries are 
ordered according to their modulus. Then set 



H = q*m*(E-e tlp S 2 -e-^S* 2 )mq 



which hence has its diagonal entries ordered according to their modulus. Setting Q = diag(g, q) 
it follows that 

-lu-loun _ ( V -1 W 1 
1 J 1 Xq*Wq 1 



Q~ l M- l SMQ 



(33) 
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Figure 2: Plot of the level spacing distribution P(s) of the unitary matrices Un (open circles) 
and , K e U N 'K e (black circles) at Lp/2-K = 0.23, A = 0.12/Vl~2 and E = 1.31 on slabs of width 
L = 52 and length N = 1000. Data for -nJJ^-Kf, with width L = 32 highlight the increasing 
agreement with the CUE distribution upon increasing the system size. Data for L = 100 shows 
P(s) at iV = 1000 only. Here and in all following figures, 100 different samples have been 
averaged over. Also, similar results have been obtained for L = 20, 32 and 42 and (p/2ir = 0.1, 
0.21, 0.25, 0.31, 0.41, 0.51, 0.61, 0.71, 0.81 and 0.91. 



where W = m*w n m, which is a self-adjoint Toeplitz matrix with first row (wq, . . . , u>£_i) where 
Wi = 2~^fc=i w n,k exp ( 2nl ^ k \ is the Fourier transform of the potential at height n. The energy- 
dependent L x L matrix \x is diagonal and real for real E. Next we want to bring the r.h.s. of 
(I33p for A = into a normal form. The characteristic equation is p 2 — fip+ 1 = for a diagonal 
complex matrix p =diag(pi, . . . , pi). We choose the first branch of the root and set 



If \pi\ < 2, then \pi\ = 1 and ^sm(pi) > and the Zth channel is elliptic. If on the other hand 
\pi\ > 2, then \pi\ > 1 and the channel is hyperbolic. Hence the projections on the hyperbolic 
and elliptic channels are given by nt = x(M > 2) and n e = xGH < 2) (here x denotes the 
usual characteristic function of an event). The matrices k and rj of Section [3] are k = pith + 
and 7] = pir e + 7Ch- If \pi\ = 1, then the eigenvalue is degenerate and k\ = 1 or Hi = — 1 and 
the channel is called parabolic; as there is only 1 eigenvector, the parabolic channel cannot be 
diagonalized and leads to a Jordan block. We exclude energies which have a parabolic block. 
As already indicated, such energies are called internal band edges |RSj . 

Next, we introduce the real diagonal matrix h{ip) = {{iih + ^ir e )(p — -)) - ^ depending on 
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Figure 3: Plot of the distribution of the eigenvalues of the unitary matrices Un (filled grey 
histogram) and n e UN'K e (solid black line) at ip/2ix = 0.23, A = 0.12/VT2 and E = 1.31 on 
slabs of width L = 52 and length iV = 1000. The dashed horizontal line at 1/27T indicates a 
normalized, uniform distribution from — 7r to 7T. 



the magnetic flux tp and then set 

jsf = ( h (<P) HvW 

V M<P)f (p + P -1 )^ + h((p)p- l -K h h(cp)±(p - p _1 )7r e + h{(f)(m h 

One verifies that this matrix is in HS(2L, C). As p, h(ip), n e and tt^ all commute, the inverse 
A/"- 1 = J*M*J is thus equal to 

jy-i = ( h (v)k(P- P~ 1 )n e + h{(p)p-K h -h(p)ir h 
\ -h{^)\{p + p- l )-n e - h^p- 1 ^ h(ip) 

Now 

K r-i{ p — i \ \r _ f P 7r h + §(p + p^W ~Up~ p" 1 ) 7 ^ 

V 1 o; v ^(p-p^k p-v+Kp+p- 1 )^ 

This matrix is the desired normal form 1Z = TZhR- e with k and 77 as defined above. Now we 
come back to (1331) and set 

T = tf- 1 Q- 1 M- 1 SMQN' . 

Then T G HS(2L,C) and 

T = Ke xv , with P = M- 1 ( ° ?V- (34) 



?W? 
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Figure 4: Plot of the distribution p(r) of the modulus of the matrix entries of n e UNiT e for the 
same parameter values as in Figure [2j For comparison the curve ( I35I) is plotted as dashed line. 
Note how the agreement between the numerical data and fl35l) gets better upon increasing the 
system size from L = 20 (L e = 13, front histogram) to L = 52 (L e = 31, back). 



As M depends on E, so does V. Note that V G hs(2L, C) and that V is nilpotent and the 
only random matrix. This is the desired normal form and (fT7|) . Now it makes sense to 
ask whether the RPP holds for the magnetic Anderson model. We do not have an analytical 
proof for the RPP, but do provide the following numerical evidence. From (I34p a realization 
of the Markov process $ n defined in (jlfip can be generated numerically (with a fixed initial 
condition $ ), giving thus a sequence of unitaries (U n , V n )i< n <N by (ITS]) . We choose N = 1000 
and thus obtain an empirical ensemble of unitaries (Un, Vn) consisting of 9000 matrices from 
100 independent disorder configurations for each of the 90 intermediate matrices produced 
during the course of the Markov process after the first 100 multiplications. As a test whether 
this ensemble is distributed as stated in the RPP, we plot the level spacing statistics P(s) for 
this ensemble (see |Meh] for a definition). Figure [2] shows the result for Un and 'kJJn'Kz for 
parameter values as stated in the figure caption. The curves obtained for Vjv and ■n e VN^ e are 
similar. Hence the results support the central property (R3) of the RPP, as does Figure [3j 
Another numerical check of (R3) is whether the matrix entries of n e UN^ e are distributed on 
the unit disc {re 10 G C | r G [0, 1], 9 G [0, 2tt)} according to the rotation invariant measure 

hlZl (i - r 2 ) Le ~ 2 rdrd6 . (35) 

7T 

The distribution of the modulus r of the matrix entries of ir e UN^e is shown in Figure HI Note, 
however, that clearly the entries are not independent. We also checked that the Hilbert-Schmidt 
norm of ir e UN^h is very small with high probability which verifies (Rl). In conclusion, we 
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consider this as sufficient evidence that the RPP holds at least approximately for the magnetic 
Anderson model. 

Now we want deduce formula for the Lyapunov spectrum from the RPP and Theorem [TJ 
Let us begin by writing out more explicit formulas for V G hs(2L,C), using the notation 

W v = h(<p)4*Wqh(<p), 



From this it follows that 



V 



n e pn e 



-n h W v % h W v ir h 



so that (V e ) 2 = and W* = W v imply 



W v n h 





-K e W v TT e 



TTeW^eWpTTe 





For the calculation of the Lyapunov spectrum using case K = C of Theorem [U we now need 
to calculate the average over V . For this purpose we use 



E(w p w q ) 



• "p+q mod L 



(36) 



This implies for any diagonal matrix d — (di, . . . , d^) 



E W*dW 



it* 

i=i / 



Thus 



ETr[n e (P*+P)n e Pn e ] = ETr(7r e W^7r e W> e 



1 

I 



\k=L h +l 



(37) 



It thus follows from Theorem [T] that, if the RPP holds, for any elliptic index Lh < p < L one 
has 



7 
'p 



4L LI 



.k=L h +l 



L-p+-)+0{\*) 



(38) 



Now replacing the definition of h(<p) and expanding in (p leads to the formula ([6]) given in the 
introduction. 
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8 Anderson model on tube geometry 



This section deals with the modifications to the prior analysis in the case of vanishing mag- 
netic field <f = 0. The transfer matrix given in ( 1321) now has real entries and S2 + 5*2 has a 
twofold degenerate spectrum (apart from the top and possibly the bottom eigenvalues). Hence 
the transfer matrices are in HS(2L,R) corresponding to the symmetry class of time-reversal 
invariant systems with even spin. All the formulas in Section [7] remain valid, but a further basis 
change is needed in order to assure that the new normal form is in HS(2L,R). All modified 
objects will carry a sombrero. 

As S2 + 5*2 i s rea l symmetric, it can be diagonalized by an orthogonal (unitary with real 
entries) which we construct first. Note that the fundamental (in our choice of sign of S2+SZ lying 
on top of the spectrum) is always non-degenerate and has always a real eigenvector fi = /x,; 
moreover, for even L, the eigenvector jhii — I'l/2 is real as well and has a non-degenerate 
eigenvalue. For other I < ^ , real normalized eigenvectors are obtained by 

(TiJl-i) = (fi,fL-i)^ (7 \ ) ■ (39) 

This is just the Fourier basis given by cosines and sinus. Now set 

m = {f u ...J L )eO{L) . 

There exists a unitary matrix a G U(L) essentially given by the r.h.s. of ( [39]) such that rh = ma 
where m is the discrete Fourier transform as defined in Section [7J If one now associates to rh as 
before a matrix M. = and carries out the associated basis change as in ([33]) . the matrices 

q and fi remain the same (hence don't carry a sombrero) and W is replaced by the real matrix 
W = a*Wa. As q and \x are unchanged, so are the diagonal matrices p, ir e , n h and g which all 
have the two- fold degeneracy (apart from the fundamentals). Also 1Z = 1Z is unchanged, but 
in ([31]) the definition of V contains W instead of W, namely V = A*VA with A = diag(a, a). 
With these changes, one has 

f = TZe xf G HS(2L,R), Pg1is(2L,R). 

Again the RPP can be checked numerically in a similar fashion as in Section [71 The results 
are of the same nature. In order to exhibit the difference between the cases with and without 
magnetic field, we plot the distribution of the level spacing of Tr e U]yTr e (Tr e VNTr e )* w 7r e [/jv(V)v)*7i" e 
for the case with and without magnetic field. For (p 7^ one has CUE statistics reflecting that 
Un and Vn are essentially independent, while for ip = the relation Vjv = Un implies COE 
statistics, both facts that can be read off Figure [5j 

Now the calculation of the Lyapunov spectrum using case K = R of Theorem [TJ we need to 
calculate ETr(n e ('P* + V) H e VTl e ). But the unitary A commutes with Il e so that the result 
is given by the formula (157jl . From this follows formula ([S]). 
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Figure 5: Plot of the level spacing P(s) of 7r e C/jv(Viv)*7re generated for the Anderson model 
at cp = (x) and tp/2n = 0.23 (•), A = 0.12/^12 and E = 1.31 on a slab of width L = 52 
and length N = 1000. The small open circles (o) show Un(Vn)*, i.e. without removing any 
hyperbolic channels. The lines denote the Wigner surmises. 



9 Ando model on tube geometry 

In this section, the symplectic universality class is dealt with. Hence we take H = H + \Hi 
with (p = and t ^ which is the Ando Hamiltonian [Ando] . Of course, the time reversal 
invariance (TRI) of if is a crucial property. Recall that for systems with odd spin (actually half 
odd spin) the time reversal operator is given by complex conjugation followed by a rotation 
in spin space by 180 degrees, which is given by the operator / = e msV = (? just as it 
was defined in Section [2j This operator acts on the spin variable only and could hence also 
be written as 1 <g> I in £ 2 (Z, C L ) <g> C 2 . The TRI property of H now reads I* HI = H which 
can readily be verified using the identity I*sl = — s. As recalled in Section [2], this means 
precisely that H has quaternion entries. This just reflects that the spin operators are linked 
to the quaternions by (q 3 , q 2 , q%) = 2zs. Similar as the objects without magnetic field carried a 
sombrero, we place here a tilde on the objects linked to the Ando model. 

The transfer matrix S associated to the Ando Hamiltonian is given in (j3J). For the analysis 
of the normal form, it is convenient to think of it as a 2L x 2L matrix with real quaternion 
entries. In fact, one may write 

~ _ f [E-(S 2 + S* 2 ) - t{S 2 - S* 2 )q 3 ](l + tq 2 )- 1 -(1 + tq 2 y \f 1 \ 

v (i+^r 1 o j \ A(i+t 2 )-^ i j ' 

' (40) 

and the TRI is then simply the fact that S G HS(2L, H) can indeed be written using the 
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quaternions with real coefficients. 

Next we want to bring S for A = into a normal form similar as in Section [71 A first 
step can be to apply the real Fourier transfer A4. already used in Section [8j As it is real, this 
conserves the TRI. The result is the analog of (|33|) : 

M-^SM = ( "™* )( , ? 1 . (41) 



i + t^r 1 o ; v A(i + t 2 )- 1 iy i 

where W is as in Section [8] and ju — fh*{E — (S 2 + S 2 ) — t(S 2 — S^)q^)m. Due to conserved 
TRI, all matrices appearing in (|4ip are in HS(2L, H). However, /J is not diagonal any more 
as it was the case without spin-orbit coupling. Indeed, it is a direct sum of (the fundamentals 
and) 2x2 quaternion blocks of the form 

E — 2 cos(r/) —2t sin (77) g 3 
2t sin (77) q 3 E — 2 cos (77) 

with a frequency 77 G ^ Z which can be calculated for each block. In (j4"Ij) this leads to 4 x 4 
quaternion blocks in HS(4, H), which are thus 8x8 matrices with complex entries. It is clear 
that diagonalizing them is a hard task. It is somewhat easier to first apply only the complex 
Fourier transform A4 = AAA*: 



M-izu - f [E-2co S (r ] )+2ttsm(r ] )q 3 }(l+tq 2 )- 1 -(l+tq 2 )* \ 

M SM - { ) 



X 



1 

A(i + t 2 )- 1 iy 1 



(42) 



Indeed, now the upper left entry and Ji — E — 2 cos(r/) + 2rf sin(^)g3 have only 2x2 blocks on 
the diagonal. These diagonals are not quaternions though (due to the imaginary coefficients 
of qs) and, in fact, the matrices on the r.h.s. are only in HS(4L, C) and not in HS(2L, H). 
On the other hand, the lowest order in A now has 4x4 complex blocks and is thus easier to 
diagonalize. It is then possible to recombine these diagonalized blocks to reinstall TRI. This 
diagonalization is carried out in detail in Appendix B. The result is the construction of a matrix 
AT e HS(2L,H) such that 

N^Q^M^SMQN = TZe xv , 

where now all factors are in HS(2L, H) and TZ is in a normal form. Therefore, one is precisely 
in the situation where the RPP for the case K = EI may hold and the Lyapunov spectrum can 
be calculated by Theorem [U In order to write out a formula for the Lyapunov exponents, one 
needs to calculate ETr(II e ('P* + V) riePIle). Because of the complicated form of Af (given in 
Appendix B) this is quite a formidable algebraic challenge, which we refrain from dealing with 
for the following reason. The matrix V and hence the coefficient in Theorem [1] depends in a 
differentiable way on t, hence for small t the coefficient is given by the one in the case K = M. 
Therefore the formula ([7]) follows. 
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10 Anderson model on a slab 



As an approximation to a <i-dimensional situation, let now consider operators on the Hilbert 
space £ 2 (Z) (gXC^)^" 1 . Here one may view (C^)^" 1 as the Hilbert space of a (d— l)-dimensional 
discrete torus and the total physical space Zx {1, . . . , iV} d_1 as a discretized and periodized 
slab. The dimension of the fiber is here L = N^ 1 . On the fiber there are cyclic shifts S2, ■ ■ ■ , 
satisfying each = 1. Given magnetic phases if 2, ■ ■ ■ , fd, the kinetic operator is now 

d 

H = St + St + A , A = ^2(e^Sj + e'^S*) . 

i=2 

The diagonal disorder is V = A J2 ni & w m with w m = J2n 2 n =1 Wn l n )( n l wnere the wjs are 
i.i.d. real and centered random variables with unit variance and n = (ni, . . . , n^). Viewing A 

also as an operator on C L , the transfer matrix now is 

5 /B-A-Auv., Its OTm al 



1 

form is found exactly in the same way as in Section [71 Let rrij for j = 2, . . . , d be the discrete 
Fourier transforms in the transverse directions and set m = m 2 <8> . . . <8> vn^. Then m is an L x L 
unitary matrix diagonalizing A. Thus also m*(E — A)m is diagonal. The ordering according to 
the modulus is again obtained by some permutation q. Then set fi = q*m*(E — A)mq. Then 
define 7r e , TCh, p, h — h((p), Af and n e using the same formulas as in Section [7J This leads again 
to a normal form with a new perturbation V. One can now carefully check that (1371) holds 
again, with h as just defined and L = iV^ -1 . Therefore Theorem [1] shows (say in the case of 
the unitary universality class) that 

Note that the number of terms in the sum over k is exactly L e and each summand is of order 
1, so that the inverse localization length is roughly equal to 7^ ~ 4^3=1 • 



Appendix A: Moments of the Haar measure on U(L) 

In this appendix, we reassemble some results, most well-known, about second and fourth mo- 
ments of the Haar measure on the unitary group U(L). The notation for the entries of a unitary 
matrix is U — (U P , q )i< P ,q<L and the average w.r.t. to the Haar measure is denote by ( . ). We 
will need formulas for the low moments of unitaries. All first and third moments vanish, and 
the only non- vanishing second and fourth moments are given by the following list (given, e.g., 
in EE): 

(Ml) (Up~ q U k j) = \ for p = k and q = I. 

(M2) (U k}P Ui ;q U m) qU n>p ) = jjzi for k = n ^ I = m and p ^ q. 
(M3) (U kjP Ui >q U mt qU n>p ) = - X(L 2_ 1} foik = m=£l = n and p ^ q. 
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(M4) (U ktP Ui !q U mi qU nt p) = L(L 1 +1) for k = n = I = m and p ^ q. 

(M5) (Uk,pUi s U mtq U rhP ) = L( - L 1 +1 - ) for either k = n^l = moik = m^l = n, both with p — q. 
(M6) {U Kp Ui s U rntq U n .p) = L(L 2 +1) for fc = n = / = m and p = q. 

For our purposes, it will be convenient to deduce the following identities. 



Lemma 2 Let A, B,C,D G Mat(L x L,C). Then 

1 



(i) ( Tr(U*AUB) ) 

(ii) ( Tr(UAUB) ) 

(iii) ( Tr(U*AUBU*CUD) ) 

(iv) ( Tr^WAUBtfCUD) ) 



(v) ( Tr(U*AUBU t CUD) ) = 



Tr(A)Tr(fi) 



= -Tr(AB* 



L 2 - 1 



Tr(A)Tr(C)Tr(BD) + Tr(AC)Tr( J B)Tr(D) 



1 



L(L 2 - 1) 
1 



Tr(AC)Tr(BD) + Tr(A)Tr( J B)Tr(C)Tr(D)J , 



L 2 - 1 



1 



Tr(A)Tr(C)Tr(BD) + Ti(AC t )Tr(BD t 

Tr(AC t )Tr(BD) + Tv(A)Ti(C)Tv(BD t ) 



L(L 2 - 1 

^(AC^TriBD*) + Tr(AC)Tr( J B)Tr(D) 



L 2 - 1 



1 



Tr(AC)Tr(5D*) + Tr(AC*)Tr( J B)Tr(D) 



L(L 2 - 1) 

Moreover, all formulas remain valid if all U's are replaced by their complex conjugates, e.g. 

{Tr^AUBtfCUD)) = (Tr(U*AUBU*CUD)} . 

Proof. Item (i) and (ii) follow directly from (Ml). Let us exemplify the proof of the others by 
treating (iii) using (M2)-(M6). First write out the trace explicitly: 

Tr(U*AUBU*CUD) = £ £ Ui^k^B^ q U^ q C m>n U n> ^D^ , 

p,q,p',q' k,l,m,n 

where all sums run from 1 to L. Now for a non-vanishing average one needs either (p', q') = (p, q) 
or (p',q') = (q,p). Hence 

(Tr(U*AUBU*CUD)) = (f&^U^U^) A k>l B p>q C m , n D q)P 

Pi^q k,l,m,n 

~t~ ^ ^ ^ ^ {Uk,pUl t qU m ^ q [J nt p) Aj t iB qq C mn Dpp 
p^q k,l,m,n 

^ ^ ^ ^ (Uk,pUl t pUm,pUn,p) A)~jBp t pC mrl Dpp . 
p k,l,m,n 
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Let us call the three summands on the r.h.s. I±, I 2 and I3. Next using (M2), (M3) and (M4) 
and the identity jt^i — l(i}-i) = lTl+T) > one nas 

II — J2 _\ ^k,kB v ^ q C m)jn D q , p — j^/jj2 _ I) Ak,mBp t qC m! kDq jP • 

Apart from the missing diagonal term p = q, this gives the summands jjzi Tr(A) Tr(C) 
Tr(BD) and l ^ j ~ 2 1 _ 1 ^ Tt(AC)Tt(B D) on the r.h.s. of (iii). Similarly, the two other summands of 
(iii) are given by I 2 again up to the missing diagonal terms. However, J 3 can be evaluated using 
(M5) and (M6). Due to the two cases in (M5) and the factor 2 in (M6), the contribution of ^3 
provides precisely both of the diagonal terms, as can be checked using once again the identity 
of fractions cited above. □ 



Appendix B: Transfer matrix of spin orbit Laplacian 

In this appendix we diagonalize 



(E — 2 cos(t?) + 2%t sin(77)g 3 )(l + tq 2 )' 1 







G HS(4,C) . 



This is one of the 4x4 blocks of (j4"2j) . Hence we may think of rj as a real number, even though 
it will be a diagonal matrix in Section [91 Note that even though is written with quaternions 
representing 2x2 blocks, it is not a quaternion matrix because there is an imaginary coefficient. 
Another important fact is that S v and S- v are similar, namely 



Q2 
q 2 



Q2 
q 2 



(43) 



Therefore, the spectra of S v and S_ v are equal and their eigenvectors linked by a simple relation. 
This then allows to build a matrix in HS(4, H) from the checker board direct sum of and 

Setting e = E — 2 cos(^) and / = 2tsin(?7), the matrix S v can be written out in a more 
explicit way 

/e-ft -et-f _■■ . \ 

1+t 2 1+t 2 1 1 \ 

et-f e+ft _. _-■ 
1+t 2 1+t 2 L L 

1+t 2 1+t 2 u u 
\ 1+t 2 1+t 2 

The characteristic polynomial is real and given by 



(44) 



/ 



p(A) = A 4 - aA 3 + b\ 2 -a\ + l = A 2 [(A + A" 1 ) 2 - a(A + A" 1 ) + 6-2] 



where 



2e 
1+t 2 



eH 2 - f 2 + 2-2t 4 
(1 + t 2 ) 2 
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Characteristic polynomials of this form always appears for real symplectic matrices and the 
stability analysis is carried out in a standard manner (as described, e.g., in |HMj ). Solving the 
quadratic equation in A + A -1 first shows 



p(\) = (A 2 -^ + A + l)(A 2 



v-X + 1) 



a a 2 



Therefore the four eigenvalues of S v are k + , — , 



k_, — where 




In these formulas the square root is taken to be the first branch so that both k± have non- 
negative imaginary part. Except for a discrete set of energies, all four eigenvalues are distinct. 
We restrict our attention to this case. One may then think of as being composed of two 
channels given by the two eigenvalue pairs k + , and k_, and the corresponding two- 
dimensional eigenspaces. If k± E S , the channel is elliptic or parabolic, otherwise hyperbolic. 

2 

Whenever ^- + 2 — b < 0, one has u + = VZ and both channels are hyperbolic but, in fact, 
not uncorrelated because k + = (TtT) -1 . Whenever ^ + 2 — b > 0, k± is elliptic if and only if 
\u±\ < 2. Resuming, the eigenvalues are generically in one of the following four constellations: 

(Gl) k + = (kI)- 1 with k + £ S 1 U R. 
(G2) k ± E SV{±1} with k + ^ «_. 

(G3) k + E §7{±1} and k_ E R/{±1} or k_ E § l /{±l} and k + E R/{±1}. 
(G4) k ± E R/{±1}. 

The cases (Gl) and (G4) are fully hyperbolic, in (G3) there is an elliptic and a hyperbolic 
channel, while (G2) is fully elliptic. 

Next let us determine eigenvectors v + ,v' + ,v_,v'_ E C 4 associated in that order to the eigen- 
values — . Due to the simple form (JUJ) of S v , the lowest two components of the 
eigenvalue equation allow to express the last two components of the eigenvectors in terms of 
the first two. Replacing this, the first two components of the eigenvalue equation become equa- 
tions in the first two components of eigenvalues only. One of them is redundant because of the 
eigenvalue property. With some care, one thus gets 



where c± E C are normalization constants to be chosen later. Using — instead of k±, the 
eigenvectors v± are constructed similarly with normalization constants c'±. As all of these are 
eigenvectors of a hermitian symplectic matrix, they have isotropy properties which we exemplify 
next (an exhaustive list is then readily made). One has {v±)*Jv± = \k±\~ 2 (S v v±)* JS r] v± = 
\k±\~ 2 (v±)*J'v±. Hence whenever k± E~ S 1 one has {v±)*Jv± = 0. Similarly {v±)*Jv T = 



/ 



2n±t - «|(et + /) \ 
2n±t 2 + n 2 ± t(et + f) 

2t — K±f — K 2 ± t 



v± = c± 



\ 



—K±e + k± J 
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whenever k±k t ^ 1, (v±)*J'v± = whenever 7€±(/t-|-) _1 7^ 1 and (v±)*J'vL = whenever 
/^(/tzp) -1 7^ 1. Note that this gives no information on {v±)*Jv± when k± G S 1 . However, the 
antisymmetry of J always guarantees that {v±)*Jv± 6 «M. Another important fact resulting 
from the reality of S v is that if «S^u± = k±v±, then also iS^t^ = 7c± v±. This may be empty if 
t>± is real (which can only happen for real k±), but produces a relation between eigenvectors 
otherwise. 

From the eigenvectors we now build a matrix M. r] g HS(4, C) such that M-'^S^M.^ is in a 
normal form. This has to be done in each of the case (Gl) to (G4) separately. 

(Gl) Here v± = vL. Furthermore the pairs (v+,v+) and (v-,vZ) span a Lagrangian plane 
(maximally hermitian isotropic for J). The normalization constants can be chosen such that 
(v + )*Jv_ = 1. Then set M v = (v + ,v'_,v^,v' + ) = (v + ,v^,v-,vZ). One can check that M. v G 
HS(4, C) and (M^^S^M^ = diag(/c+,7q:, (k+) _1 , (T^)" 1 ) G HS(2, H). Below we will use that 
also the complex conjugate A4 V G HS(4, C) diagonalizes S v . 

(G2) Here the eigenvalues are e llf+ 1 e~ lif+ 1 e lip ~ 1 e~ l{f ~ with eigenvectors v+,v+, v-,vZ. Neither 
of these eigenvectors are isotropic. Thus we introduce the real vectors w± = 9?e(t>±) and 
w'± = $tm(v±). Hence all the vectors w±,w± are isotropic (because {w±)*Jw± is both real 
and imaginary). Furthermore {w±)*Jw T = and {w±)*Jw' T = 0. By adequate choice of 
the normalization constants one can achieve (w±)* J~w'± = 0. All this assures that M. v = 
(w + ,W-,w' + ,w'_) G HS(4,R). Furthermore, one can check 



(M v )~ SjjMjj 



f cos(ip + ) — sin((/9 + ) \ 

cos(y9_) — sin((^_ 

sin(<£> + ) cos(y? + ) 

V Sin(<y2_) COs(</2_) J 



G HS(4,C) . (45) 



(G3) This is a combination of (G2) above and (G4) treated next. 

(G4) Here all eigenvectors are real and thus isotropic, and, moreover, one has {y±)*JvL = 0, 
{y±)*Jv T = and {y' ± )*Jv' T = 0. Hence {v+)*Jv- and {y' + )*Jv'_ cannot vanish (otherwise 
there would be a three-dimensional isotropic subspace) and can by choice of the normalization 
constants both be made equal to 1. Then J\4 V = (v + , v'_, v' + , V-) G HS(4, R) and one can verify 
{M-r^^S^M.^ = diag(/c+, («_) _1 , G HS(2, C). Note that the complex conjugate of 

the latter diagonal matrix is also in HS(2, C). 

In each of the four cases let us denote T> v G HS(4, C) the matrix such that 



(M n y S v M n = v v , (M n y S V M V = V 



v ■ 



Note that the complex conjugate equation differs from the first one only in case (Gl). Further 
by®, ~ 

(M- ri y 1 S- v M- v = V„, = ( q * M M v . 



q 2/ 

We apply the checker board sypmplectic sum to the matrices S v and S- v with quaternion 
entries (but with complex coefficients), and consider also the 4x4 matrix A4 V as a 2 x 2 matrix 
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of 2 x 2 blocks so that also Ai v (BAi- v is well-defined. Then the above translates to 



and each factor is in HS(8, C). Now consider similar as in Section [8] the matrix A = diag(a, a) 
with a 4 x 4 matrix a = 2~v (~?\), having diagonal 2x2 blocks. Then A*S V ®S_ V A G HS(4, H) 
and 

(A*^T v ®M- v By l (A*S ri ®S- v A) (A*M~ n ®M- v B) = {B- l V~ n @V n B) , 
with B G HS(8, C) chosen such that all four factors are in HS(4, H): 



B = diag(6,&) 



/l 0\ 

10 

10 

\0 1/ 



In case (G2), B^V^QV^B is given by the r.h.s. of ( 1451) . albeit each entry multiplied by go- The 
cases (G3) and (G4) are similar, but in case (Gl) one needs to use that already V v G HS(2,H). 
As to A*J^i v ®Ai- v B, it is in HS(8,C) because each factor is, and then it is tedious but 
straightforward to show that it actually has quaternion entries. 



Finally resuming, let us set ftf n = A*M V ®M- V B, then 

(Sf v )- 1 (A*S v ®S^ n A) = B^W^V^B G HS(4, 
From this one readily builds the matrix M needed in Section |9j 
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